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ABSTRACT 

High-quality observations of B and V light curves obtained at Las Campanas Ob- 
servatory for local Type la Supernovae (SNe la) show clear evidence that SNe la with 
the same brightness decline or stretch may have systematic and independent deviations 
at times t < 5 days before and at times t > 30 days after maximum light. This sug- 
gests the existence of two independent secondary parameters which control the shape of 
SN la light curves in addition to the primary light curve parameter, stretch s or Amis. 
The secondary parameters may reflect two independent physical effects caused by vari- 
ations in the initial carbon-to-oxygen (C/O) profile in the progenitor and the initial 
central density p c in a carbon-oxygen white dwarf exploding as a SN la. Theoretical 
light curves of delayed detonation SN la models with varying progenitor masses on the 
main sequence, varying accretion rates, and varying primordial metallicity reproduce 
two morphologically different and independent types of variations in observed visual 
light curves. These calculations predict small variations of ~ 0.05 mag in the abso- 
lute brightness of SNe la which are correlated with the variations of progenitor mass 
on the main sequence Myis which changes the C/O profile, and p c which depends on 
the accretion rate. Such variations in real supernovae will induce systematic errors in 
SN la calibration at high redshifts. A physically motivated three-parameter, s, C/O, p c , 
template for SNe la light curves might take these variations into account. Comparison 
between the theoretical predictions and the observational results agree qualitatively; 
however, the observations show variations between the B and V light curves that are 
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not expected from the modelling and may indicate limitations in the details of the 
theoretical models. 

Subject headings: Supernovae: general 



1. Introduction 

SNe la are thought to be thermonuclear explosions of massive carbon-oxygen white dwarfs 
(CO-WD) in binary stellar systems. These supernovae are important tools ("standard candles") 
of modern physics and cosmology. Use of SNe la as standard candles has provided th e first direct 



evidence of the accelerating e xpansion of the universe and the existence of dark energy (jRiess et al 



19981 : IPerlmutter et al.lll999T ). 



Maximum luminosity varies among SNe la and is not a constant. Essentially all supernova- 
based cosmology studies use methods of calibration of SNe la as standardizable candles (removing 
luminosity variations) which rely on empirical relations between the intrinsic brightness at maxi- 
mum light and other observable characteristics of SNe la such as the shape of the light curve, or 



the rate of brightne ss decline after maximum light - the brightness-decline relation (|Phillips!ll993l ; 



Phillips et al.lll999l ). The present accuracy of calibration, ~ 10%, has been sufficient for discovering 
the Dark Energy, but it must be improved to perhaps ~ 1 — 2% in order to study properties of 
Dark Energy quantitatively. This is a formidable task which requires increasing the accuracy of 
SN la light curve observations, accounting for effects of dust absorption, and so on. It also requires 
improving the calibration procedure itself, which involves two important interconnected issues. 

(1) If SN la light curves form a one-dimensional family characterized by a single parameter 
such as Ami5, then the observed spread of individual SNe la around the average brightness-decline 
relation can be attributed to random statistical error, and the accuracy of cosmological measure- 
ments can be be increased by simply increasing the number of observations of individual SNe la. 
On the other hand, if light curves are characterized, in addition to Amis, by some yet unknown 
independent "secondary" parameters, then improving the calibration is impossible without taking 
the dependence of the light curve on secondary parameters into account. So far, it is unknown if 
secondary parameters exist. Clear evidence for at least two independent secondary parameters will 
be provided in this paper. 

(2) Implicit in all SN la calibration procedures is a fundamental assumption that nearby 
and distant (cosmological) SNe la behave identically and that empirical brightness-decline or 
brightness-stretch relations established for local SNe la can be used for high-redshift cosmolog- 
ical supernovae as well. Obviously, empirical studies of nearby and distant SNe la alone can not 
confirm or reject the existence of variations of brightness-decline relations with redshift. This 
requires independent accurate measurements of the intrinsic brightness of SNe la. Secondary pa- 
rameters may hold a key to this difficult problem. If we understand the physical mechanisms and 
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the relationship of secondary parameters to initial conditions of SN la explosions, e.g., metallicity 
and/or the parameters of binary progenitors, we may gain some insights into systematic changes 
in SN la lightcurves with cosmological time. 

The Carnegie Supernova Project has recently obtained a highly uniform set of SNe la light- 
curves with an accuracy of a few hundredths of a magnitude both for individual SNe l a and in 



terms of variations between different objects (jContreras et al.ll2009l ; iFolatelli et al.l 12009). These 



new data provide clear evidence for the existence of secondary variations in SNe la light curves 
which are independent of the primary Phillips relation, and thus, are evidence for the existence of 
independent secondary parameters. The CSP data allow us to begin addressing issues (1) and (2) 
outlined above. 

SN la models predict a weak dependence of the early light curve on the carbon-to-oxygen 
(C/O) ratio of the progenitor. They also predict a weak dependence of the late-time light curve 
on the initial central density, p c , of the exploding WD. A combined analysis of new observational 
data and theoretical predictions leads us to suggest the existence of two independent secondary 
parameters that, in addition to Amis ° r stretch, control the intrinsic brightness of SNe la. 

The paper is organized as follows: § Ogives a short review of previous work on this subject; §[3] 
describes new observational data used in this paper; § 0] briefly summarizes the current theoretical 
understanding of SNe la and presents theoretical calculations of SNe la light curves; § [5] analyzes 
the observational light curves which give the evidence for the existence of secondary parameters, 
compares observational and theoretical light curves, and discusses the theoretical interpretation and 
the mechanisms by which secondary parameters arise in SNe la; results of the paper are summarized 
and the discussion of the implications of the results to the calibration of SNe la is presented in § [6j 



Previous Work 



Over the last half decade, a number of observational and theoretical studies have sought to 
uncover secondary parameters. Much of this effort has been in attempts to fin d direct correlations 



between physical effects and peak lum i nosity. Exampl e s include metallicity (IWang et al 



Hofiich et al 



1998; 



Piro k Bildsten 



Howell et al 



Timmes et al 



2008 



2001 



2001 



20031: lEllis etaDbood : [Gallagher et al 



Kasen et al 



2003 



2008 



Chamula k et al.ll2008r). asymmetri e s of the explosion (IWang et al 



Kasen et al 



2004; 



Ropke et al 



2006 



Hofiich et al 



Hoflichl 120061 ) . ag e of the progenitor 



1997a 



Howell et al.l 



2006; iKasen et al 



2009 
1997b| 
2009; 



Wang Wheeler 20od\ . central densi ty and C/O ratio ( Hofiich et al. 1998 . 200ol; IPommguez et al 



Man nucci et alj|2006l ; ISullivan et al 
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20061 ). neutron-rich isotope to 56 Ni rat io ( Mazzali k Podsiadlowskil 120061 ) . and the opacity of the 



overlying material ( Mazzali et al. 2001 : Kasen k Woosleyl 2007 ). Our approach here is different. 
Rather than focusing on a single physical effect on a specific observable or stage, we make use 
of detailed stellar evolution m odels that were calculated all the wa y through white dwarf forma- 
tion, accretion, and explosion (IDommguez et al.ll200ll : iHofhchl 120061 ) and focus on the effects that 
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variations in the progenitor and accretion rate have on both the early and late parts of the light 
curve but using parametrizations when the underlying physics is uncertain. Both approaches are 
valid, and we must take care not to draw conclusions beyond the range of validity of our 1-D 
approach and rely on consistency checks using observations. At the present time, this approach 
makes sense since 1-D models do a reasonable job of reproducing observations. While current 3-D 
models are able to reproduce some of the observations, no 3-D model to date has begun with a 
confirguration that was the result of detailed stellar evolution, nor is it understood what physi- 
cal variations in 3-D models reproduce the the tight brightness-decline relation. With or without 
constant mixing, 1-D models can reproduce the brightnes s-decline relation and its narrow width 



(Hoflich et al. 



1996 



Umeda et al.lll999l ; iHoflich et al.l 12002 ). a varying amount of mi xing produces 



an 'anti-correlation' and/or a huge spread comp arable to the entire range of SNe la (jHoflich et al 



1996 



Pinto Eastman 



2000 



Kasen et al.l 120091 ). 



T hough details depend on the pre-condit ioning of the WD (|Hoflich &: Steinl 120021 : iLivne et al 



2005 



Kasen et al. 



2009 



Zingale et al.1 120091 ). 3D models of deflagrations and delayed detona- 



tions predict strong mixing of the cent ral region during the deflagration phase (jKhokhlov 



Gamezo et al. 2003; 



2004; 



Gerardv et al 



^ein ecke et al.ll2002l ) in conflict with obse rvations of late time spectra ( Hoflich et al 



2007]) and remnants (IFesen et al.l 120071). We note t hat recollapsing m odels 



200i 



0; 



(jBravo &: Garcia-Sen z 2006; Baron et al.ll2008l ; iBravo Garcia-Senzll2009l ; IBravo et al.ll2009l ) avoid 



central mixing, but produce only 56 Ni in the center. In these models the center has expanded sum 
ciently prior to carbon ignition that the electron capture rates hav e dropped and so burning occurs 



without neutronization. The gravitational confined detonation (IPlewa et al.1 12004 : iJordan et al 



2008: lMeakin et al.ll2009l ) model will p roduce neutrqnized material under some ignition conditions 
and avoid its production under others (IJordan et alJ 12009). 



3. Observations 

The light curves of SNe are often stitched together from observations carried out on a variety of 
telescopes at a variety of sites. This has the advantage of filling in gaps in the light curves. However, 
there is a distinct disadvantage. The spectral energy distributions of SNe are considerably different 
than those of normal stars, and the spectra of SNe change on time scales of days. While some 
spectral features are easily associated with singly and doubly ionized metals such as silicon and iron, 
other absorption features are actually blends of many lines. The effective bandpasses of filters vary 
from camera to camera. The net result is that photometry of SNe carried out on different telescopes 
often exhibits systematic offsets, sometimes amounting to 0.2 mag. From synthetic photometry of 
spectra of normal stars and spectra of SNe at different times with resp ect to maximum light 



we can compute "S-c orrections" which largely resolve these differences (jStritzinger et al 



2002 



Krisciunas et al.ll2003l ). However, unless we have good sequences of spectra for all of our SNe, 
it is not possible to devise error- free S-corrections. Certainly, we wish to attribute variations in 
light curve morphology to the SNe themselves, not to some conspiracy of the telescopes, sites, and 
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cameras. 



The Carnegie Su pernova Project (C SP), which began operation in September 2004, seeks to 
address this problem (jHamuv et al.l 12006). The CSP endeavors to observe Type la and Type II-P 
SNe in the filters of the Sloan Digital Sky Survey (u'g'r'i'), the standard Johnson B and V filters, 
plus the near-infrared bands Y, J, and H. Roughly 50 SNe are being followed each year, for five 
years. Almost all of the optical photometry of nearby SNe is being obtained with the Swope 1-m 
telescope at Las Campanas Observatory (LCO). For higher redshift SNe some BVI data are derived 
from images with the Wide Field CCD camera on the 2.5-m DuPont telescope at LCO, and a small 
amount of BVR data is from the 6.5-m Clay Telescope (Ma gellan #2). Based on h igh-quality CSP 
data the RMS deviations of 0.025 mag have been achieved (jContreras et al.l l2009). 



We use 18 SNe la observed in 2004, 2005, and 2006 as part of the CSP. All these had well 
sampled light curves. Many were observed a week or more prior to maximum light. Many were 
observed 60 or more days after maximum light. For the purposes of this paper we restrict ourselves 
to the B— and, in particular, V— band photometry because, based on theoretical models, V is 
expected to be least sensitive to additional variables. 



The BVRI templates of SNe la are provided in lPrieto et al.l ( 20061 ). Sets of BVRI templates 
are characterized by the standard dec li ne ra te parameter Amis (2?) which serves as a morphological 
label for the set. While iPrieto et al.1 (120061 ) provide BVRI templates on l y BV were used in this 
analysis. Prior to fitting the light curves with the templates of lPrieto et al.l (|2006l ) we first estimated 
the time of B-band maximum and subtracted the redshift- and time- dependent K-c o rrecti ons from 
the photometry by interpolating the B- and F-band corrections of IHamuv et al.l (|1993l ). If the 
subsequently determined time of maximum light was more than 0.5 days different than the value 
adopted for the calculation of the K-corrections, then they were recalculated. It was then checked 
that the newly determined time of maximum light was statistically consistent and thus there was 
no need to iterate further. 



The range of Amis(-B) for the IPrieto et al.l (|2006l ) te mplates is 0. 8 3 — 1 .93. For example, for 



a given object we chose B- and F-band templates from IPrieto et al.l ((2006), stretched them by 
the time dilation factor (1+z), then shifted them by small increments over a range of dates and 
over a range of magnitudes to minimize the total x 2 °f the fit in each filter. Then we tried the 
other templates including all bands over a range of Anii^(B) to determine which templates give the 
lowest x 2 °f an - I n this wav we determined the observed maximum magnitudes and the times of 
maximum light. For each of our SNe we obtained functions that fit the B- and l/-band light curves 
from roughly -5 until 25 days after the time of maximum light. This time interval for fitting the 
data was motivated by the theoretical models which predicted that, during this peri od, the visual 



LCs s hould be least effected by variations in th e central density and the progenitors (jrlonich et al 



1998 



Dommguez &i Hoflich 



2000 



Hoflicbl l2006). 



For a comparison of the fits to the photometry, we then subtracted off the derived maximum 
magnitudes in B and V, subtracted off the times of maximum light in the two filters, and divided the 
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"time-since-maximum" by (1+z) to give "rest frame days since maximum" . The same normalization 
is applied to the photometric data. 



4. Models of SNIa 



4.1. Explosion 

Current observations of SNe la favor a delayed detonation (DD) scenario of a SN la explosion in 
which the explosion begins as a subsonic deflagration which later tur ns into a supers onic detonation 



by the process of a deflagration-to-detonation transition or DDT (IKhokhlovl 1 1 99 ll ) . The ensuing 



detonation incinerates the entire WD. In one-dimensional models the deflagration speed, S, and the 
transition density, pt r , at which the DDT occurs are free parameters. The value of p± r determines 
the fraction of a WD that will burn to nuclear statistical equilibrium and produce 56 Ni. Therefore, 
ptr is a critical parameter that controls the brightness of a SNe la. DD models with pt r — 0.5 — 2.5 x 
10 7 g cm -3 reproduce the observed range of SNe la luminosities, correct stratification of chemical 
elements in SNe la envelopes, and the cor relation between maximum brightness and width of 



2001 



Hdflichetal 



Marion et a 



(Fesen et al 



200 



3L 



2002 



Hoflich 



2006 



the l i ght curve consisten t with observations (Hoflich 1995: Hoflic h et al. 1996: Mazza li et al 



Quimby et al.1 120071 : iMazzali et al J 120071 ; iKasen et al 



1998, 



2009 



2007 



DD models appea r to be in agreement with observations of SN la remnants 



Badenes et al J 12008 ). 



The initial central density of a WD, p c , its metallicity, and its C/O ratio also influence the 
production of 56 Ni and maximum brightness of SNe la during the explosion. Note that in our 
models, the free parameters are the m ain sequence mass of the progenitor, Mmsj the accretion 



rate, and the primordial metallicity Z (IDommguez et al 



2001 



Hoflich! l2006h . The fact that the 



brightness of a SN la is controlled, for the most part, by a single, primary parameter pt r , and 
that there exists an evolutionary bottleneck associated with the limiting Chandrasekhar mass of a 
WD provides a plausible explanation for a near one-dimensional sequence of SN la explosions. A 
crucial point is that p c , primordial metallicity, and C/O ratio influence certain characteristics of the 
explosion not associated with the primary brightness-decline relation. In particular, the C/O ratio 
has an influence on the expansion velocity of SNe la. The dependence is caused by variations in 
nuclear binding energy of the CO fuel. The larger the C/O ratio, the smaller the binding energy, and 
the faster the SN la envelope expands. This is a secondary effect since the bulk of the kinetic energy 
is determined by a much larger difference in binding energies of C/O and products of explosive 
burning (Fe-peak elements). On the other hand, p c influences the distribution and amount of 56 Ni 
in the innermost parts of a SN la. SN la models based on explosions of Chandrasekhar-mass WDs 
predict a hole in the 56 Ni distribution near the center which is filled instead with highly neutronized 
isotopes of Fe-group elements. The hole is caused by electron captures and neutronization of matter 
at high densities. The size of the hole increases with p c . These two effects influence the formation 
of both early and late portions of a light curve. As noted above the nickel hole is generally absent in 
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all 3-D models, either because the core expands prior to ignition, or because in delayed detonation 
models mixing occurs during the deflagration phase, but the stratified structure is restored by the 
detonation. Because the flame is very topologically complex during the deflagration phase, the 
detonation will not produce a central neutronized hole. 

The initial C/O ratio and profile and p c in a supernova are a result of stellar and binary 
evolution. The C/O profile is produced during the central Helium burning and thin shell burning 
during the stellar evolution and the accretion to M h- The central He burning is initially dominated 
by carbon production via the 3a reaction in the convective core. When the He mass fraction 
becomes depleted, 12 C(a, 7) 16 mainly controls He-burning and most of the 16 is produced during 
the late phases of central He-burning. Note that the final abundances depend on a combination of 
the 12 C(a,7) 16 rate and chemical mixing which determines the duration of the phase of depleted 
He-core burning. Stellar evolution models which produce low C/O ratios of ~ 0.25 — 0.4 are 
in agreement with observational constraints, like the a mount of oxygen found in the inner zone 



of pulsating white dwarfs or the age of open c lusters (jDominguez et al. 



2001 



Metcalfe et al 



2001 



1999 



Dommguez et al 



Straniero et al.1 120031 ) . The C/O ratio in the burning shell is greater, 



~ 1, because the shell helium source has lower density and higher temperature compared to helium 
burning in the core. The size of the convective core depends mainly on the progenitor mass Mms ° n 
the main sequence and, to so me extent, on the primordial metallicity Z, n amely the iron abundance, 



which dominates the opacity (jDommguez et al. 



2001 



Hoflich et alJl200d ) The dependence of M c 



on Mms is nonlinear, with the mass of the He convective core changing with Mms slowly for 
low Mms an d more rapidly when Mms approaches its maximum value. Approximately, ]VI core is 

0. 3.0.4.0.7M for M M s = 1.5,5,7M , respectively (jDommguez et alJl200lh . The core's C/O ratio 
is smallest for small MmSj increases rapidly when Mms increases from 1.5 to 3M and remains 
relatively constant for larger Mms- These two effects combine to make the mean C/O ratio in an 
exploding WD decrease with Mms slowly for low mass stars Mms — 1-5 — 5M and much more 
rapidly for higher mass stars with Mms — 5 — 7M . In turn, this makes the explosion energy of a 
WD more sensitive to Mms when Mms is large. 

The mass of the WD core is sensitive to the mass of the progenitor on the main sequence and is 
rather insensitive to the primordial metallicity. To first order, the C/O ratio and thus, the explosion 
energy decreases with Mms • This tendency is generally true for Mms larger than 3M independent 
of assumptions about mixing and nuclear reaction rates. For lower mass progenitors, the total mean 
C/O ratio varies little with Mms- M core decreases with Mms, but this effect is almost compensated 
by a decrease of the local C/O ratio which decreases with Mms- Both this effect and the high C/O 
ratio in shell burning can be understood by nuclear physics and chemical mixing. The local C/O 
ratio depends mostly on the competition between triple-a, 3 4 He — > 12 C, and a-capture on 12 C, 

1. e. 12 C(a,7) 16 0. As a three-body reaction, triple-a dominates at high a concentrations whereas 
12 C(a,7) 16 dominates when 4 He is depleted (and at high temperature). During the early phases 
of helium core burning mostly 12 C is produced but, eventually it is depleted by 12 C (a, 7) 16 O. For 
burning in small convective cores, even moderate chemical mixing keeps 4 He at a certain level and 
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prolongs the phase of 4 He burning under depleted conditions, reducing carbon. In burning in thin 
radiative shells, the temperatures are higher and the burning time-scales are shorter, therefore less 
carbon is depleted. 

The central density p c at which the WD ignites is controlled by the competition between 
adiabatic compression caused by accumulation of mass at the surface and energy losses in the 
center of the WD. As a result, p c depends on a configuration of the binary system, that is on the 
rate of accretion onto the WD from the stellar companion. 



The physics of the ignition process is multi-D in nature. With the exception of lGarcia-Senz fc Woosley 



()1995l ) who find ignition occurs in rising plumes, all mult i-D simulations to date show ignition at or 
near the center due to the downward motion of plumes (jHoflich fc Steinll2002l : IZingale et al.l l2009) 
in the simmering phase. Though relevant for the preconditioning of the explosion, variations in the 
final stage of runaway are of the order of hours (|Hoflich fc Steinl |2002| ) and unlikely to effect the 
central density. As a result, p c depends on the configuration of the binary system, the evolutionary 
state of the stellar companion, and the resulting accretion rate. 

In what follows, we characterize the initial conditions by p c , and the progenitor characterized 
by the initial metallicity Z and Mms- We use initial distributions of C and O in the core as 
predicted by the evolutionary calculations of a main sequence star with appropriate Mms- Both 
during stellar shell burning and He burning during the accretion C/O ~ 1. Once the stellar and 
WD evolution has been calculated to the onset of the explosion, the explosion is calculated using 
one-dimensional DD models taking into account the progenitor evolution, the hydrodynamics of 
the explosion, detailed nuclear networks with 213 isotopes, the radiation t ransport, detailed atomic 
mod els, and 7-ray transport. Details of the actual models are described in lDommguez et al.l (|200ll ) 
andlHoflichl feoodi ). 



Dommguez et al.1 (120011) described the effects of varying the C(a,7) O rate from the high 



value (jCaughlan et al.l Il985l ) used in our calculations to the lower value of ICaughlan &: Fowler 
(|1988l ) . The effect on the final comp ositions is large, but we use the value that has been shown to 



agree with observational constraints (IDominguez et a l. 1999; Dommguez et al.ll200ll ; lMetcalfe et al 



200 ll ) Note, however, that the final C/O ratio depend s on both the rate an d the chemical mixing 
scheme adopted in the stellar models (see for example I Straniero et al.ll2003l ). 



4.2. Theoretical light curves 



We calculated B and V light curves of a series of DD models with fixed pt r = 2.3 x 10 7 g cm -3 
which produce explosions with spectral and light curve characteristics of normal bright SNe la. 
Our fiducial model has Pr = 2 x 10 9 g cm -3 , primordial metallicity Z equal to the solar value Z@ 
(jAnders k, Grevessdll989l ). and Mms = 5M . The abundances in the WD are a result of the stellar 
evolution of a main sequence star with meta llicity, Z^. When sc aling Z, for elements up to Si, we 
adopted the [0 /Fe] abundance suggested by lArgast et al.l (|2000l ) which implies smaller variations 
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with redshift for elements up to Si compared with Fe by a factor of three. This is done because 
22 Ne affects the explosive nucleosynthesis whereas Fe determines the opacity, and therefore, the size 
of the convective He-burning core (for a given mass) and, to some extent, the -B-band magnitudes. 
Note that the effect of primordial metallicity Z on the explosive nucleosynthesis is dominated by the 
22 Ne abundance and not by the iron abundance because it is the reduction of the proton/nucleon 
ratio, Y e , which changes the explosive equilibrium abundances. The light curves presented below 
illustrate the effects that variations of central density and progenitor mass, Mms, have on light 
curves of SNe la. 



4-2.1. Influence of progenitor mass Mms an d metallicity 

Fig. [1] presents B and V light curves and B—V of four DD models with fixed p c = 2x 10 9 g cm -3 
and varying Mms an d primordial metallicity (upper panels). The lower panels present a differential 
comparison of light curves normalized to maximum light. The difference is defined as dM(t) = 
M(t) — Mf(t), where My is the magnitude of the fiducial DD model. The figure shows that 
variations in the progenitor, i.e. the main-sequence mass Mms an d metallicity, strongly change the 
rise to maximum light. These variations are caused by variations in expansion velocity in models 
with various Mms and hence with various C/O ratios. The expansion velocity decreases when Mms 
increases and the overall C/O ratio decreases. At the same time, the chemical and density structure 
of the outer parts of the SN la envelope is similar for all DD models. Therefore, variations in the 
formation of the early light curve are mostly controlled by the rate at which the outer layers expand 
and become transparent. The faster the expansion rate, the faster the photosphere recedes, and 
the faster the light curve rises towards maximum (see the upper left panel of Fig. [1]). The light 
curve of the Mms = 7M Q model in Fig. [1] rises notably slower than the fiducial light curve while 
the light curve of Mms = 1-5Mq model rises somewhat faster. This is reflected as the early time 
negative and positive alM in the left panels of Fig. [1] The effect is more pronounced for higher Mms 
due to the greater sensitivity of M core on Mms near the high end of Mms interval. In particular, a 
strong secondary extremum develops when M core extends to layers which only undergo incomplete 
silicon burning. As a consequence, we expect that differential effects are most pronounced in stellar 
populations with a mix of young and old stars, since these will contain a distribution of Mms- 
We notice also that variations of the progenitor lead to small variations on a 10 percent level in 
B — V and its evolution with time (upper right panel). Changing the primordial Z will increase 
the primordial iron abundance in the outer layers, and, at the expense of 56 Ni, more 54 Fe will be 
produced from 22 Ne. Primordial metallicity plays a minor role for variations in V , but B — V 
becomes bluer with decreasing primordial metallicity . This direct 'phq tospheric' effect does not 



change the visual LCs but the B and U V light curves (Hoflich et al.lll99 8). The evolution in B — V 



is similar to and has been discussed in iKrisciunas et al.l (|2003l ) 



If we were to "observe" these four supernovae, match their light curves by applying the stretch 
correction, and then compare the residual differences we would get the result shown in the lower 
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right panels of Fig. [TJ We show dM(t) after matching the light curves around maximum and 
determining the stretch over the range up to ~ 15 days past maximum light. The end result is 
essentially identical light curves near and past maximum with the exception of a strong deviation in 
the rise time for the light curve with the highest Mms = 7M Q and a small deviation of the same light 
curve approximately 30 days past maximum. The pre-maximum deviation of the Mms = 1-5M 
light curve is smaller but still visible at the level of ~ 0.1 mag. 

There has long been a s uspicion that the metallicity of the progenitor should be associated 
with the luminosity at peak (jWang et al.l l200ll ; lHamuv et al.l Il995l ; iBranch et al.l Il99a ) . Recently 
attempts have been made to measure directly the averag e metallicity in the en v ironment of the 



SN u sing either line indices or ratios (jHamuv et al.l |2000| ; iGallagher et al.l 120051: iGallagher et al 



2008) or by measuring the spectral energy distribution (SED) of the galaxy ([Howell et al.ll2009i ). 
These studies haved failed to confirm the expected shift in peak luminosity with metallicity. Nu- 
clear physics predicts that primarily due to increased amounts of the neutron rich 22 Ne there 
should be a direct correlation between the amount of 56 Ni p roduced and the progenitor metallic- 



ity (iTimmes et al.l 120031 ) . Recent results ([Howell et al 



200C 



Neill et al 



expected dependence of metallicity on -W. where,, balWherlTTl £53) report a metallicity 



2010) find a weaker than 



depen dence in the Hubble residual, but not the peak luminosity. From the results of lArgast et al 



(|2000l ) the variation of elements below Si (including Ne) varies less than Fe by about a factor of 
three, thus the sensitivity of the peak brightness is smaller than would be found by just scaling all 
elements to the iron abundance. Theoretical attempts to find correlations between metallicity and 
peak luminosity (or light curve shape) have not followed the detailed stellar evolution through to 
explosion, but have rather just altered the nucle osynthetic yields p ost-explosion and calculated the 
light curve, scaling just on the iron abundance (jKasen et al.l l2009). 



4-2.2. Influence of central density p c 

Fig. [2] shows the same type of comparison as Fig. Q] but for a series of models with fixed 
Mms = 5Mq and solar metallicity, and with varying accretion rate which leads to a central density 
p c = 1.5 x 10 9 , 2 x 10 9 , and 6 x 10 9 g cm -3 which, in our models, corresponds to late-time accretion 
rates of 1 x 10~ 7 — 2 x 10 -6 M yr" 1 . Fig. [2] shows small differences between the light curves prior 
to maximum light. They are virtually unnoticeable in the upper panels of Fig. [2] but can be clearly 
seen on the differential plots (lower left panel). These differences arise due to small variations in 
the binding energy of WD models with different p c . When p c increases, so does the binding energy, 
and this translates into a somewhat smaller expansion velocity. As a result, light curves with 
larger /smaller p c rise slower /faster and this results in a negative/positive pre-maximum differential 
dM. 

We can also see that variations in p c have a significant effect on the behavior of light curves 
which begin to show up at ~ 20 — 25 days after maximum light. Variations in p c lead to a noticeable 
shift of the late time light curve with respect to the absolute magnitude at maximum light. The 
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light curves shift up or down when p c decreases or increases, respectively. This effect is related 
to the existence of the central hole in the distribution of 56 Ni in Chandrasekhar-mass models of 
SNe la. Due to increasing electron capture with density, the nuclear statistical equilibrium shifts 
away from 56 Ni to stable isotopes of the iron group when the central density of a WD increases. 
Near maximum and shortly past maximum light the envelope of a SN la is rather opaque and 7-rays 
emitted near the center are trapped and do not contribute to the formation of the light curve. The 
light curve around maximum light is controlled by the distribution of 56 Ni in the outer parts of 
the supernova. As time goes on, the envelope expands and the distribution of 56 Ni near the center 
begins to influence the formation of the light curve. Larger central density means a larger hole and 
less 56 Ni. As a result, the light curve of a SN la with higher p c becomes shifted down with respect 
to maximum. With decreasing p c the hole is smaller which means more 56 Ni and the resulting shift 
of the light curve is positive. Note that B — V is very similar but at late times differs from the 
'general blue shift' produced by variations in progenitors. The lower right panels of Fig. [2] again 
show the differentials when a stretch correction has been applied. Small pre-maximum differences 
in the light curves have virtually disappeared. However, shifts in the late light curves caused by 
variations of 56 Ni near the center remain very pronounced. 

Variations in some of the characteristics of light curves associated with variations in p c and 
Mms are summarized in Table 1. Table 1 shows that brightness of SNe la in both My and Mb 
increases monotonically from —19.25 to —19.11 and —19.32 to —19.18 mag, respectively, when 
Mms decreases from 7 to 1.5M when other parameters, pt r , Z, and p c , are kept constant. These 
variations are of the order ~ 0.1 mag. Variations of My and Mb with p c are smaller, ±0.03 mag, 
because p c mostly affects the electron-capture in the center of the SN la which hardly contributes 
to the SN la luminosity at maximum. The effect of metallicity Z can be seen by comparing the 
fiducial model (first model) of Table 1 with Z = 0.02 and the fifth model with Z = 0.002. Z affects 
both Mb and My and has a pronounced influence on the B — V colors of SNe la. 

Variations in p c and Mms have pronounced secondary differential effects on stretch-matched 
light curves as illustrated in Fig. [3l (1) Mms influences the rise time of the light curves prior to 
maximum light. Larger Mms leads to slower rise and vice versa. The upper right plot of Fig. [3] 
shows secondary variations of the absolute visual magnitude of SNe la, My, and differentials in My 
as a function of Mms f° r times 25, 42, and 55 days past maximum light. (2) On the other hand, p c 
influences the light curve ~ 30 days after maximum and later. Increases or decreases in p c cause 
the later portion of the light curve to shift down or up with respect to maximum, respectively. The 
upper left plot of Fig. [3] shows My, and differentials in My as a function of p c at 20 and 40 days past 
maximum light. Secondary variations in B — V for both series of models are given as a function of 
Mms ° r Pc in the lower left plot of Fig. [3j All the plots illustrate the point that there are noticeable 
secondary variations of color and absolute visual magnitude of stretch-matched supernovae which 
are associated with variations of secondary parameters. 

Finally, the lower right plot of Fig. [3] summarizes the relation of a predicted relative variation 
in absolute brightness of SNe la (the quantity which cannot be determined from observations unless 
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the absolute brightness of SNe la has been measured using an independent method of calibration) 
and relative secondary variations (differentials) of stretch-matched light curves of SNe la which 
may be directly accessible to observations. This plot shows that secondary variations in My of 
stretch-matched supernovae may reach 0.2 mag. 



5. Analysis of Observations 

The above theoretical considerations guide us in our differential analysis of observed light 
curves of SNe la. 

Below we present a differential comparison of V and B light curves of several SNe la obtained 
by the CSP survey (Figs. H]-[8|). The basic LC properties are given in Tabled These objects 
were selected because they were discovered well before maximum, have good time coverage for 
up to ~ 60 days past maximum light, and very small observational errors. SN 2005na serves as 
the reference or fiducial model for the comparison. Stretch corrections for all SNe la are listed 
in Table 1. Left panels on each figure show a comparison of light curves with stretch corrections 
from Table 1. Right panels show the comparison of SNe la after we added an additional stretch 
to their light curves in order to match the s-factor of the fiducial model, SN 2005n a . Th at is, 



rather than just making the Amis stretch correction using the formulae of iJha et al.l (|2006l ). the 
curves were stretched such that the luminosity was brought into coincidence with SN 2005na. The 
additional stretch is rather small since all these objects are in the range of normal bright SNe la. 
This procedure helps reduce primary differences in due to the brightness decline relation (see §H]). 

SN 2005al vs. SN 2005na (Fig. 0]). These supernovae are well within the normal bright range. Early 
portions of their light curves, less than 25 days after maximum light, are very similar. This is also 
indicated by the very small difference in s-factors of these supernovae (Table [2]) . Exact matching 
gives only a marginal improvement in dispersion in M — M re f. Fig. 0] shows that M — M re f is less 
than a few hundredths of a magnitude until ~ 25 days after maximum. At later times the V light 
curves are shifted with respect to each other by ~ 0.2 mag. There is no discernible systematic shift 
in B. This behavior in V might be explained by variations in p c and the size of 56 Ni hole in the 
center of the SN la. 

SN 2004ef vs. SN 2005na (Fig. [5]). These two supernovae are different in their pre-maximum 
behavior and show an extremum between 20 and 40 days in V, both characteristics of variations in 
the progenitor. There might be a small systematic shift, ~ 0.05 mag, in the late portions of the V 
light curves as well, although it is less pronounced compared to that of the previous pair SN 2005al 
and SN 2005na. Differences between SN 2004ef and SN 2005na may indicate both variations in 
progenitor masses and metallicity, and rather similar p c in these two events. 

SN 2005ki vs. SN 2005na (Fig. [6]). Both supernovae are very similar early on with a very small 
shift in the late time light curves of less than 0.1 mag. They likely have similar values of p c , Mms> 
and metallicity. 
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SN 2005el vs. SN 2005na (Fig. [7]) These two supernovae show differences in the early light curves, 
a secondary extremum, and a significant shift in the late time light curves as well. Note that 
the stretches of SN 2005el and SN 2004ef are very similar, making it not very likely that the huge 
M(SN 2005el)— M re f can be attributed to differences in s. In terms of SNe la models, the difference 
hints at variations in progenitor mass M^s an d p c - 

SN 2005am vs. SN 2005el (Fig. [8]). Finally, we present a differential comparison of two super- 
novae with a large difference in s- values. SN 2005am, is a steep decliner with Amis = 1-61 
whereas SN 2005el is a normal SNe la with Amis = 1.37. After stretching (Fig. [8]), the values of 
M(SN 2004am) and M re f are surprisingly similar. Within the framework of models, this suggests 
similar progenitor masses and p c . 

Contrary to theoretical calculations which predict similar behavior of M — M re f in both B 
and V, observations show, depending on the example, variations in the B and V morphology. At 
the moment the reason is not clear. Several effects may play a role: a) errors in s and the time of 
maximum t max , b) differences in the theoretical and observational filter functions B and V, and c) 
limitations of the explosion models and progenitors. 

First, we investigate the stability of differential comparisons of light curves. Matching the time 
of maximum light t max of two supernovae may introduce an error in the differential comparison. 
Another source of the differential error may be small variations in Amis or s - 

Fig. [9] shows the effects of a relative shift of t max and a variation of Amis of supernovae 
SN 2004am and SN 2004el. These two supernovae represent a pair with different rise times and a 
corresponding deviation at ~ 30 days after maximum. As discussed in § 0] this may be attributed 
to variations in the progenitor masses of these two objects. The differentials between supernovae 
of various progenitors are more sensitive to uncertainties in the observables than those with the 
p c -signature because the former show variations at both early and late times. 

Maxima were shifted by ±1 day and 2 days, and s by ±0.1. The figure clearly shows that the 
morphology of M — M rc f in V is stable, namely the rise, a rather flat part around maximum, and 
a dip at day 30. The variation of the differentials over the dip is less than 0.05 mag. For B, we 
still see that the main characteristics of M — M re f in B are stable, namely extrema at ~ 15 and 40 
days, but they are larger than those in V in size, e.g., the first dip varies by almost 0.3 mag. 

The difference in the shapes of B and V differentials appear to be real including the differences 
in B between observed and predicted morphology. In part, the size of the first £>-dip may be caused 
by fringing, i.e. shifts of LCs with several maxima and minima, but hardly goes away. 

One possible effect is a wavelength offset in the B filters. Another group of effects includes 
limitations inherent to the models: a) Metallicity variations for stars with the same progenitor 
show differences of less than 0. 05 mag (Fig. [T|), a value consistent with detailed spectral analyses 



(jLentz et al.ll200ll ; lHonichll2006l ). However, uncertainties in our progenitor evolution during central 



He-burning may underestimate the variations in progenitors, b) In B at this level of accuracy, we 
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may be seeing limitations inherent in spherical models. After maximum light, B is very sensitive to 
temperature variations because it is formed by the Wien tail of the source function, S oc e ~ hu / kT ^ 
where S is the source function in the radiative transfer equation. For example, small variations of 
56 Ni mixing will increase the temperature and, at the same time, the blocking in B. Indeed, off- 
cente r DDT models show spectral changes i n B which are compatible to the size of the deviation we 



find (Hoflich et al 



2006; 



Kasen et al.ll2009l ). However, strong rotational mixing is rather unlikely 



because the impact on the brightness decline relation and other observational constraints discussed 
in § [H As seen above, uncertainties in the determination of t ma x and s have moderate influence on 
the differential in V, but may amplify differences in B. 



6. Discussion and Conclusions 

We analyzed a set of high quality uniform V and B light curves of Type la super novae ob- 
tained by the Carnegie Supernova Project. These data provide clear evidence for the existence of 
secondary variations in SN la light curves which are independent of the primary Phillips relation 
and, thus, for the existence of at least two independent secondary parameters. Comparison of the 
data with a series of non-LTE light curve calculations of delayed detonation explosions indicates 
that these secondary parameters may be physically related to variations of central density, p c of 
a Chandrasekhar-mass white dwarf exploding as a SN la and to variations of the main sequence 
mass, MmSj of the primary stellar companion in a progenitor binary stellar system. 

It is generally accepted that the total amount of 56 Ni produced in the explosion is the pri- 
mary parameter which controls the absolute brightness and the rate of post-maximum decline. 
Production of nickel depends mainly on the transition density pt r at which subsonic burning turns 
into a supersonic detonation. However, if p± r were the only parameter describing the explosion we 
would have a perfect one-dimensional sequence of SN la events and a one-dimensional family of 
light curves. Some previous attempts to find additional parame ters were done by comparing the 



risetimes of the nearby SNe la sa mple to that of the dist ant one. iRiess et al.l (| 19991 ) claimed to see 



a significant difference, however, iGoldhaber et al.l (|200ll ) did not find the effect to be statistically 
significant. 

The high quality light curves of SNe la obtained by the CSP clearly show that one parameter 
is insufficient to characterize the light curve. Light curves of several CSP supernovae shown in 
§ [5] illustrate the fact that two supernovae with identical behavior at maximum light and the early 
portion of the post-maximum light curve, — 30 days past maximum, may behave differently prior 
to maximum as well as at late times (> 30 days past maximum light). In particular, SNe la that 
are similar near maximum may have different pre-maximum rise and may have the late portion of 
the V light curve (> 30 past maximum light) shifted up or down by different amounts with respect 
to maximum. Examples presented in § [5] show that these differences may reach 0.2 — 0.3 mag. 

Light curve modeling (§ UJ predicts that two supernovae with identical behavior of light curves 
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near maximum may exhibit small differences in the rise time to maximum and small shifts of the 
light curves with respect to maximum luminosity at times > 30 days after maximum. These vari- 
ations may be traced to variations in initial central density, p c , and the WD progenitor mass on 
the main sequence, Mms- Instead of Mms, one may consider the carbon-oxygen ratio C/O in the 
pre-exploding WD as a second parameter. As far as light curves of SNe la are concerned, these 
parameters, p c and Mms (or C/O), may be treated as "independent" secondary parameters char- 
acterizing the light curve in addition to its main characterization via Amis or stretch. Theoretical 
predictions seem to be qualitatively consistent with two distinct morphological types of deviations 
shown by observations in V (§ ED - The B light curves also show distinct morphological types of 
deviations but the morphology appears to be different from that of the V light curves. On the 
contrary, light curve modeling predicts similar morphological behavior in V and B. We note that 
most of the energy flux is coming in the V band and thus the calculations of the V light curves 
should be more reliable. The reason for the discrepancy in B is not clear and requires further 
investigation. 

The existence of independent secondary parameters in SNe la has several important implica- 
tions. Obviously, any calibration of supernovae using a one-parameter set of light curve templates 
should lead to systematic errors in template matching and to calibration errors. The results of this 
paper indicate that a set of light-curve templates should form at least a three-parameter family. 
Our analysis suggests that a physically motivated set of templates may be constructed by using a 
primary SN la template parameterized by Amis ° r stretch with (1) an additional correction of the 
slope of the pre-maximum light curve and (2) an additional offset of the late V light curve with 
respect to maximum. Physically, the primary parameter Amis should be thought of as a parameter 
which reflects variations in the amount and distribution of 56 Ni in the outer parts of the SNe la 
envelope responsible for the formation of the early post-maximum part of the light curve. The rise 
time correction might reflect the variation in the C/O ratio, and the offset of the light curve at 
later times might reflect variations in initial central density p c and the amount and distribution of 
56 Ni in the central parts of the envelope. 

However, producing a multi-parameter set of light-curve templates is not sufficient for improv- 
ing the calibration procedure. The crucial point is that variations of initial conditions responsible 
for secondary parameters must also cause, according to theoretical predictions, small variations in 
absolute brightness of SNe la. The calibration must take the dependence of the absolute brightness 
on secondary parameters into account. By matching the proposed mult i- parameter templates with 
high-quality observations it should be possible to account for secondary variations in the intrinsic 
brightness and to reduce the calibration errors. High quality data are not sufficient at this time for 
carrying out such a program systematically. However, this work clearly indicates the potential of 
high quality uniform sets of observations for studying secondary variations in SN la light curves and 
for providing important theoretical clues about physical mechanisms of such variations. One may 
hope that with improvements in SN la theory, verified with observations of nearby supernovae, it 
would be possible to predict variations of absolute brightness of SNe la as a function of secondary 
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parameters. This, in turn, may provide an opportunity to predict evolutionary effects in SNe la 
calibration related to systematic changes in initial conditions with redshift. 

Finally, we have to address the promise and limitations of this study. Our quantitative theo- 
retical predictions depend on our exact treatment of the explosion mechanism and stellar evolution. 
Thus, variations in the treatment of convection, the 12 C(q,7) 16 rate, rotation, and 3-D explosive 
effects could be important at some level. Nevertheless we have tested the trends that we predict on 
a small number of SNe la and they should be tested on a large and homogeneous set of SNe la. Our 
preliminary results are promisin g. Light curves for a la rge number of SN ar e/will be available b y 
proje cts such as the E SSENCE dESSENCdboogh . CFH (jSNLSlbood ). NSNF (|Aldering et alJbood ). 
PTF |Rau et al-lhoolih . and LSST (|LSSTil200gh . and we will publish an extensive comparison. The 
properties of the components, i.e. the shape in M — M re f, may be based on theoretical models 
and optimized using large data sets. However, even if our results are confirmed fully, systematic 
theoretical studies including 3-D effects will be essential to go further. As discussed in § [H SN la 
physics is intrinsically 3-D and thus those effects must be taken into account. However, the last 
decade or so of theoretical work indicates that the effect should not be dominant since it appears 
from the the observations that 3-D effects like rotation of t he WD and the position of initial igni 



tion are reduced by the effects of the deflagration and DDT ( Khokhlov 199a: iKhokhlov et al 



Niemever 



199' 



Livne et al 



al: IrS 



einecke et al 



2002; 



Gamezo et al. 



20051 : iRopke Hillebrandt 



2005 



2003; 



Gamezo et al.l 120051: iPlewa et al 



1997 



2004 



e.g.). For a review see iHoflichl (|2006l ). Finally, we 
return to the implications of asymmetry for the use of SNe la for cosmology. A 10% asymmetry of 
the photosphere would not cause systematic difficulties for us ing SNe la as distance indicators at 
the current level of accuracy of about 20% (|Wang et al.ll2003l ) though it would require that most 
of the dispersion has to be attributed to asymmetry. We note that if such effects are present, then 
SNe la are even more homogeneous than they seem from current dispersions in peak brightness. 
This level of asymmetry w ould, however, cause a directional dependence of the luminosity of order 
~ 0.1 mag (jHoflichl Il99ll ) and a corresponding, but smaller, dispersion in the brightness-decline 
relation of SNe la. This dispersion depends on the viewing angle dependence of the luminosity 
variation and, thus, the nature of the asymmetry. The angle dependence of the luminosity due to 
the viewing angle 9 of a single SN la will not, in general, vary as the line of sight to the equator 
as cosO. A mor e stringent limit co mes from observations of individual supernovae. The first broad 
band survey by lWang et al.l (119961 ) established that SNe la have very low polarization at a level of 
P ~ 0.2% whereas core-collapse SNe are generally more highly polarized P ~ lquality polarimetry 
measurements indicating that SNe la are more highly polariz ed before maximum, and a few weeks 
past maximum the polarization disappears (jWang et al.ll2003l ). SN 2004dt showed that SNe la have 
spectral feature dependent polarizati on, implying different chemical species have different geometry 
(jWang et al.ll2006l ; iPatat et al.1 12009). Since the continuum polarization at maximum light is ob- 



served to be P < 0.2% which, for scattering dominated atmospheres, translat e s into a directional 



dependence of the flux at the 0.05 % lev el (Hoflich 



2001 



Wang fc Wheeler 



2004 



Wang et alJ 12006 : 



1991 



Wang et alJll9 97b. 



Fesen et al 



2007; 



Patat et al 



2003 



Howell et al 



2009), one expects 



that while 3-D effects are important for understanding the explosion mechanism, their observational 
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effects are relatively small. However, as summarized in I Wang Sz Wheeled (|2008l ) asymmetries lead 
to a dispersion in the color terms which can be magnified significantly when extinction corrections 
are applied. Other independent factors will contribute to the error and dispersion such that the 
interaction within the progenitor system or the primordial metallicity require early time spectra, 
and/or a combination of optical and IR data. 
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Fig. 1. — B (blue) and V (red) and B — V are given for a number of delayed detonation models 
with the same explosion parameters ( Dommguez et al.ll200ll ) and p c = 2 x 10 9 g cm -3 (see §0]) but 
varying progenitor mass Mms between 1.5 and 7.0-M© and metallicities Z between 0.002 and 0.02 
(solar). Models are referenced by the pair of numbers [Mms,-^] for [1.5,0.02] (dashed), [7.0,0.02] 
(dotted) and [5.0,0.002] (dash-dotted). The reference model has M MS = 5M , and Z = 0.02 
(solid). The B (blue) and V (red) magnitudes and the color index B — V are given in the upper 
right and left panel respectively. The Amis f° r both V and B light curves are close to within 0.03 
mag but they are not identical. The lower panels show the B and V differentials without and with 
stretch-correction to the Sf of the reference model on the left and right, respectively (see Table 1). 
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Fig. 2. — Same as in Fig. [I] but compares models with fixed Mms = 5M and Z = 0.02 and with 
varying central densities p c /(10 9 g cm" 3 ) of 1.5 (dashed), 2.0 (solid), and 6 (dotted). The central 
density is due to variations in the accretion rate (see Fig. 6 of iHoflichl 120061 ) . 
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Fig. 3. — Basic observables in V for models with varying central densities p c (crosses) and main 
sequence progenitor masses Mms (open circles) relative to the reference model. All quantities 



are n ormalized to the same fiducial s factor using the brightness decline relation (jHonich et al 



2002). The change in absolute brightness at maximum light is AM*(V), the differential brightness 
is M — M re f at times +t after maximum light for the p c -series (crosses, upper left) and the I s * 
(fa +7 days), 2 nd (« +35 days), and 3 rd + 10 days for the Mjvis-series (open circles, upper right). 
The corresponding colors B — V (lower left) and M — M re f (AM* (V)) (lower right) are shown. 
Note that the sign of the residual between two SNe la is arbitrary depending on the choice of a 
reference object. 
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SN2005al vs. reference SN2005na 
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Fig. 4. — Comparison of SN 2005al to SN 2005na chosen as a reference supernova. Left column 
- comparison of un-normalized supernovae. Right column - comparison of supernovae normalized 
to the maximum brightness in V and with the s-factor of SN 2005al adjusted to be equal to the 
s-factors of SN 2005na The original s-factors of all supernovae are listed in Table O Upper row - 
B and V light curves. Middle row - differential in V. Lower row - differential in B. 
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SN2004ef vs. reference SN2005na 
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Fig. 5. — Same as Fig. [5] but comparison of SN 2004ef to SN 2005na. 
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SN2005ki vs. reference SN2005na 
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Fig. 6. — Same as Fig. [5] but comparison of SN 2005ki to SN 2005na. 



-28- 



SN2005el vs. reference SN2005na 
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Fig. 7 — Same as Fig. H but for SN 2005el. Note the S-shape in M_M ref with a 'spread' of about 
0.15 mag. 
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SN2005arn vs. reference SN2005el 
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Fig. 8. — Same as Fig. [7] but comparison of SN 2005am and SN 2005el. Here, we used an offset 
in stretch s and in ts of 0.01 and 0.5 days compared to Table 2, respectively. The morphology of 
M — M re f does not change, but the graph is 'tilted' by about 0.1 mag. 
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Fig. 9. — Influence of uncertainties in Amis an d the time of maximum t raax on the differential 
comparison of SN2004ef and SN2005na (Fig. [SJ. Panels are marked by the assumed variations in 
Amis and t max , (5(Ami5),5(t max )). Time is in days. The functional form of the differential in 
V appears to be stable including the early rise, the extended plateau and the extremum at about 
3 weeks past maximum. The functional form of the differential in B also appears to be stable 
although quantitatively the differential is more sensitive to the variations. 
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Table 1. Properties of Calculated SNe la 



Parameter Model A Model B Model C Model D Model E Model F Model G 



M MS 


5.0 


1.5 


3.0 


7.0 


5.0 


5.0 


5.0 


Z 


0.02 


0.02 


0.02 


0.02 


0.002 


0.02 


0.02 


Pc 


2.0 


2.0 


2.0 


2.0 


2.0 


1.5 


6.0 


My 


-19.21 


-19.25 


-19.22 


-19.11 


-19.15 


-19.19 


-19.26 


tv 


18.24 


18.12 


18.19 


19.5 


18.52 


18.24 


18.24 


B-V 


-0.02 


-0.06 


-0.03 


-0.01 


-0.13 


0.03 


0.01 


M B 


-19.23 


-19.32 


-19.27 


-19.18 


-19.28 


-19.22 


-19.26 


S/Sf 


1.00 


0.96 


0.97 


1.02 


0.98 


1.0 


0.99 



Note. — Characteristics of theoretical light curves of DD models with various progenitor 
mass Mms/Mq, metallicity Z, and central density p c in units of 10 9 g cm -3 . Model A is the 
fiducial model with M M s = 5.OM , Z = 0.02 and p c = 2.0 x 10 9 g cm" 3 with Am 15 = 1.25. 
Listed for all models are the absolute maximum brightness M B / V , time of V maximum ty 
(in days), and a correction to a stretch parameter, s/sf, required for making the s-factor of 
the model equal to that of the fiducial model, Sf. 
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Table 2. Properties of observed SNe la 



OAT 

SN 


Amis 


s° 


» / c 

tB/mB 


i 1 c 

tv/mv 




5V[r/s) a 


2004ef 


1.47 




264.96 


264.54 


0.09 / 0.05 


0.00 / -0.03 




1.33 


0.88 


16.92 


17.06 


0.03 / -0.07 


0.14 / 0.04 


2005al 


1.19 




429.47 


430.96 


-0.01 / -0.02 


-0.01 / -0.01 




1.24 


0.92 


15.08 


15.08 


0.01 / -0.03 


0.16 / 0.13 


2005am 


1.56 




437.10 


437.53 


0.12 / -0.01 


0.02 / -0.07 




1.61 


0.75 


13.76 


13.84 


0.14 / -0.09 


0.29 / -0.06 


2005el 


1.36 




646.86 


647.51 


0.14 / 0.07 


0.04 / 0.00 




1.37 


0.86 


15.24 


15.22 


0.13 / 0.02 


0.24 / 0.08 


2005ki 


1.44 




705.98 


706.20 


0.12 / 0.03 


0.03 / -0.01 




1.41 


0.85 


15.69 


15.65 


0.10 / -0.04 


0.24 / 0.07 


2005na 


1.19 




740.32 


741.79 


0.00 / 0.00 


0.00 / 0.00 




1.19 


0.95 


16.26 


16.25 


0.00 / 0.00 


0.00 / 0.00 



a Top values of Amis are the values for a family of BV templates for 
which the B-band template gives the best fit to all £?-band data. Bottom 
values are derived from the "early" subset of B-data. This subset includes 
data which extends to the inflection point, some time 20 — 25 days after 
maximum light, where the second derivative of changes sign. 

b s-factor is derived from the bottom value of Amis (for the early sub- 
set of B-data; s ee note (1)) using relation between s and Amis given in 



Jha et al.1 (|200J) 



c Time (top number) and apparent maximum brightness (bottom num- 
ber) at the truncated Julian Date is "JD - 2,453,000". 

d Average difference (in mag) between B and V light curves relative to 
SN 2005el. Pairs of numbers are differences for raw (r) and stretched (s) 
data. Top pair is determined using the early subset of data. Bottom pair 
is determined using all data. See note (a). 



